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MOTION ESTIMATION PROCESS AND SYSTEM 
USING BLOCK-MATCHING AND INTEGRAL PROJECTION 

BACKGROUND OF THE INVENTIQNJ 

The present invention relates generally to methods and apparatus for 
motion estimation for video image processing, and in particular, is directed to 
improved methods and apparatus for determining motion vectors between 
frames of video images using a sparse search block-matching motion 
estimation technique and integral projection data. 

Advancements in digital technology have produced a number of digital 

J 

video applications. Digital video is cun-ently used in digital and high definition 
TV, videoconferencing, computer imaging, and high-quality video tape 
recorders. Uncompressed digital video signals constitute a huge amount of 
data and therefore require a large amount of bandwidth and memory to store 
and transmit. Many digital video systems, therefore, reduce the amount of 
- digital video data by employing data compression techniques that are 
optimized for particular applications. Digital compression devices are 
commonly refen-ed to as "encoders"; devices that perform decompression are 
referred to as "decoders". Devices that perform both encoding and decoding 
are referred to as "codecs". 

In the interest of standardizing methods for motion picture video 
compression, the Motion Picture Experts Group (MPEG) issued a number of 
standards for digital video processing. MPEG-1 addresses digital audio and 
video coding and is commonly used by video devices needing Intermediate 



data rates. MPEG-2 is used with devices using higher data rates, such as 
direct broadcast satellite systems. • 

Motion picture video sequences consist of a series of still pictures or 
"frames" that are sequentially displayed to provide the illusion of continuous 
motion. Each frame may be described as a two-dimensional array of picture 
elements, or "pixels". Each pixel describes a particular point in the picture in 
terms of brightness and hue. Pixel information can be represented in digital 
form, or encoded, and transmitted digitally. 

One way to compress video data is to take advantage of the 
redundancy between neighboring frames of a video sequence. Since 
neighboring frames tend to contain similar information, describing the 
difference between frames typically requires less data than describing the 
new frame. If there is no motion between frames, for example, coding the 
difference (zero) requires less data than receding the entire frame. 

Motion estimation is the process of estimating the displacement ■ 

' . , ■ ■ ■ 

between neighboring frames. Displacement is described as the motion ; ; 
vectors that give the best match between a specified region in the current 
frame and the corresponding displaced region in a previous or subsequent 
reference frame. The difference between the specified region in the current 
frame and the corresponding displaced region in the reference frame is 
referred to as "residue". 

In general, there are two known types of motion estimation methods 
used to estimate the motion vectors: pixel-recursive algorithms and block- 



matching algorithms. Pixel-recursive techniques predict the displacement of 
each pixel iteratively from corresponding pixels in neighboring frames. Block- 
matching algorithms, on the other hand, estimate the displacement between 
frames on a block-by-block basis and choose vectors that minimize the 
difference. 

In conventional block-matching processes, the current image to be 
encoded is divided into equal-sized blocks of pixel information. In MPEG 
video compression standards, the pixels are grouped into "macroblocks" 
consisting of a 16x16 sample array of luminance samples together with one 
8x8 block of samples for each of the two chrominance components. The 
16x16 array of luminance samples further comprises four 8x8 blocks that are 
typically used as input blocks to the compression models. 

Fig. 1 illustrates one iteration of a conventional block-matching 
process. Current frame 120 is shown divided into blocks. Each block can 
be any size, however, in an MPEG device, for example, current frame 120 
would typically be divided into 16x16-sized macroblocks. To code current " 
frame 120, each block in current frame 120 is coded in terms of its difference 
from a block in a previous frame 1 1 0 or upcoming frame 1 30. In each 
iteration of a block-matching process, cun-ent block 100 is compared with 
similar-sized "candidate" blocks within search range 1 15 of preceding frame 
1 1 0 or search range 1 35 of upcoming frame 1 30. The candidate block of the 
preceding or upcoming frame that is determined to have the smallest 
difference with respect to current block 100 is selected as the reference block. 
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shown in Fig. 1 as reference block 150. The motion vectors and residues 
between reference block 150 and current block 100 are computed and coded. 
Current frame 120 can be restored during decompression using the coding for 
each block of reference frame 110 as well as motion vectors and residues for 
5 each block of current frame 1 20. 

Difference between blocks may be calculated using any one of several 
known criterion, however, most methods generally minimize error or maximize 
correlation. Because most correlation techniques are computationally 
intensive, error-calculating methods are more commonly used. Examples of 
'1 0 error-calculating measures include mean square error (MSE), mean absolute 

distortion (MAD), and sum of absolute distortions (SAD). These criteria are 
described in Joan L. Mitchell et al., MPEG Video Compression Standard. ■ ? 
International Thomson Publishing (1997), pp. 284-86. SAD is a commonly : 
used matching criterion. . , _ 

15 . 

SAD is defined as: . , ■ ^ : 

SAD{i,J)= 5^ r(x.y)-s(x+/,y+y) | 

where block size is MxN, r(x.y) is the current block and s(x+i,y+j) is the 
candidate block within a search area 1 1 5 in the reference frame. The motion 
vector is the value (i,j) that results in the minimum value for SAD(i,j). 



A block-matching algorithm that compares the current block to every 
candidate block within the search range is called a "full search". In general. 
larger search areas generally produce a more accurate displacement vector, 
however, the computational complexity of a full search is proportional to the 
size of the search range and is too slow for some applications. A full search 
block-matching algorithm applied on a macroblock of size 16x16 pixels over a 
search range of + N pixels with one pixel accuracy, for example, requires 
(2xN+1)2 block comparisons. For N= 16. 1089 16x16 block comparisons are 
required. Because each block comparison requires 16x16, or 256, 
calculations, this method is computationally intensive and operationally very 
slow. Techniques that simply reduce the size of the search area, however, 
njn a greater risk of failing to find the optimal matching block. 

As a result, there has been much emphasis on producing fast 
algorithms for finding the matching block within a wide search range. Several 
of these techniques are described in Mitchell et al., pp. 301-11 . Most fast 
search techniques gain speed by computing the displacement only for a 
sparse sampling of the full search area. The 2-D logarithmic search, for 
example, reduces the number of computations by computing the MSE for 
successive blocks moving in the direction of minimum distortion. In a 
conjugate direction search, the algorithm searches in a horizontal direction 
until a minimum distortion is found. Then, proceeding from that point, the 
algorithm searches in a vertical direction until a minimum is found. . Both of 



these methods are faster than a full search but frequently fail to locate the 
optimal matching block. 

Another method for reducing the amount of computation in a full 
search is to calculate the displacement between blocks using integral 
projection data rather than directly using spatial domain pixel information. An 
integral projection of pixel information is a one-dimensional array of sums of 
image pixel values along a horizontal or vertical direction. Using two 1-D 
horizontal and vertical projection arrays rather than the 2-dimensional array of 
pixel information in a block-matching algorithm significantly reduces the 
number of computations of each block-matching. This technique is described 
in a paper by I.H. Lee and R.H. Park entitled "A Fast Block Matching 
Algorithm Using Integral Projections." Proc. Tencon '87 Conf., 1987, pp. 590- 
594. 

Fast motion estimation techniques are particularly useful when 
converting from one digital video format to another. Digital video is stored in 
encoded, compressed form. When converting from one format to another 
using conventional devices, the digital video must first be decompressed and 
decoded to its original pixel form and then subsequently encoded and 
compressed for storage or transmission in the new format. Conversion 
techniques requiring that digital video be fully decoded are very time- 
consuming. 

The present invention provides improved methods and apparatus for 
the motion estimation process by performing a fast search that minimizes the 



number of block comparisons while maintaining the quality of the motion 
vector. In addition, the present invention provides methods and apparatus 
for motion estimation using the fast search process of the present invention 
and integral projection to further minimize the number of computational 
operations. The present invention further provides methods and apparatus 
for fast motion estimation using integral projection that allow digital video data 
conversion from one format to a second format without full decoding to pixel 
data thereby greatly reducing the time required for data format conversion. 



J 



SUMMARY OF THE INVENTION 

Methods for obtaining a motion vector between two frames of video 
image data consistent with this invention comprise the steps of: (a) 
selecting a two-dimensional current block of a first frame, the current block 
comprising an N x M array of pixel information, where N and M are positive 
integers greater than 1; (b) determining a two-dimensional search area of a 
second frame based on the current block and a predetermined search range; 
(c) determining a first set of candidate blocks by selecting at least one 
candidate block in each strip of pixel information based on the cun-ent block 
and a predetermineci difference criterion, each candidate block comprising an 
N X M array of pixel information; (d) determining a second set of candidate, 
blocks by selecting at least one candidate block in each of at least one .; 
secondary search areas based on the current block and the predetermined 
difference criterion, each candidate block comprising an N x M array of pixel 
infomnation and each of the at least one secondary search areas based on 
the first set of candidate blocks; (e) obtaining a reference block from the ;| " . 
second set of candidate blocks based on the predetermined difference 
criterion; and (f) determining a motion vector representing the distance 
between the current block and the reference block. 

A motion estimation system for obtaining a motion vector between two 
frames of video image data comprises means for selecting a two-dimensional 
current block of a first frame, the current block comprising an N x M array of 
pixel information, where N and M are positive integers greater than 1 ; means 
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for determining a two-dimensional search area of a second frame based on 
the current block and a predetermined search range: means for determining a 
first set of candidate blocks by selecting at least one candidate block in each 
strip of pixel information based on the current block and a predetermined 
difference criterion, each candidate block comprising an N x M array of pixel 
information: means for determining a second set of candidate blocks by 
selecting a candidate block in each of at least one secondary search areas 
based on the current block and the predetermined difference criterion, each 
candidate block comprising an N x M array of pixel information and the at 
least one secondary search areas based on the first set of candidate blocks; 
means for obtaining a reference block from the second set of candidate 
blocks based on the predetermined difference criterion: and means for 
determining a motion vector representing the distance between the current 
block and the reference block. 

A computer program product consistent with the present invention 
comprises a computer-usable medium having computer-readable code 
embodied therein for obtaining a motion vector between two frames of video " 
image data, the computer-usable medium comprising a component 
configured to select a two-dimensional current block of a first frame, the 
current block comprising an N x M array of pixel information, where N and M 
are positive integers greater than 1: a component configured to determine a 
two-dimensional search area of a second frame based on the current block 
and a predetermined search range; a component configured to determine a 



first set of candidate blocks by selecting at least one candidate block in each 
strip of pixel information based on the current block and a predetermined 
difference criterion, each candidate block comprising an N x M array of pixel 
information; a component configured to determine a second set of candidate 
blocks by selecting at least one candidate block in each of at least one 
secondary search areas based on the current block and the predetermined 
difference criterion, each candidate block comprising an N x M array of pixel 
information and each of the at least one secondary search areas based on 
the first set of candidate blocks; a component configured to obtain a 
reference block from the second set of candidate blocks based on the 
predetermined difference criterion; and a component configured to determine 
a motion vector representing the distance between the current block and the 
reference block. 

According to another aspect of the invention, a method for converting a 
frame of digital video data from a first format to a second format, comprising 
the steps of obtaining digital coefficient information representing pixel 
information of a first frame in a first format; and determining a motion vector 
representing the difference between a first frame in a second format and a 
second frame in the second format based on the digital coefficient 
information, wherein the first frame in the second format corresponds to the 
first frame in the first format. 

Furthermore, a system for converting a frame of digital video data from 
a first format to a second format comprises means for obtaining digital 
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coefficient information representing pixel information of a first frame in a first 
format; and means for determining a motion vector representing the 
difference between a first frame in a second format and a second frame in the 
second format based on the digital coefficient information, wherein the first 
frame in the second format corresponds to the first frame in the first format. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The accompanying drawings, which are incorporated in and constitute 
a part of this specification, illustrate the Invention and, together with the 
description, serve to' explain the principles of the invention. 

Fig. 1 is a diagram illustrating a prior art block-matching technique; : 

Fig. 2 is a diagram illustrating a step of selecting a search area for a 
block-matching search consistent with the present invention; 

Figs. 3A and 3B are a flow diagram illustrating a method for performing 
a block-matching search consistent with the present invention; 

Fig. 4 is a diagram illustrating a step of performing an x-directional .. ; 
search of a block-matching search consistent with the present invention; v: ' 

Fig. 5 is a diagram illustrating the step of performing a y-directional 
search of a block-matching search consistent with the present invention; 

Fig. 6A is a diagram showing a method of calculating a vertical integral 
projection for an 8x8 an-ay of pixel information; 

Fig. 6B is a diagram showing a method of calculating a horizontal 
integral projection for an 8x8 array of pixel information; 
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Fig. 7 is a diagram showing a candidate block horizontally straddling 
two 8x8 blocks consistent with.the present invention; 

Fig, 8 is a diagram showing a candidate block vertically straddling two 
8x8 blocks consistent with the present invention; 

Fig. 9 is a diagram showing one step of a block-matching process 
consistent with the present invention; 

Fig. 1 0A Is a diagram showing a method of calculating vertical integral 
projections for an 8x8 array of pixel information; 

Fig. 1 0B is a diagram showing a method of calculating horizontal 
integral projections for an 8x8 array of pixel information; 

Fig. 1 1 is a diagram of a system consistent with the present invention. 

Fig. 12 is a diagram of a processor consistent with the present 
invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Reference will now be made in detail to preferred implementations ^^i- 
consistent with the present invention, an example of which is illustrated in the 
accompanying drawings. 
A. Sparse Search Process 

Generally, motion estimation techniques take advantage of 
redundancy between successive frames of a video sequence. Often, the 
neighboring frames contain very similar information although slightly 
displaced. Motion estimation techniques compress the amount of data 
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needed to represent a digital video sequence by encoding one frame in terms 
of its difference from a neighboring frame rather than encoding each frame in 
its entirety. When the video sequence is replayed, the decoder reconstructs 
the current frame using the reference frame as well as the motion vectors and 
the residues. 

In a method consistent with the present invention, a motion vector is 
estimated for each block of a current frame with respect to a reference frame. 
Fig. 2 shows current frame 220 and reference frame 210. Current frame 220 
Is divided into a plurality of blocks 225, each comprising M x N pixels. M and 
N may represent any integer values greater than one, however, in the 
following example, M=N=8. To represent current frame 220, each block 225 
in current frame 220 is represented as a motion vector indicating 
displacement from a similar size block in reference frame 21 0 and is 
represented as residues between the block in the current frame and the block 
in the reference frame. 

Figures 3A and 3B contain a flow chart illustrating a method consistent 
with the present invention for obtaining a motion vector for one block of a : : 
current frame. In particular, Figs. 3A and SB describe the method of 
estimating a motion vector for one block. To begin, the process selects 
current block 240 of current frame 220 as shown in Fig. 2 (Step 302). To 
determine a motion vector using a block-matching algorithm, current block 
240 is compared with some or all of the candidate blocks in a con-esponding 
search area of reference frame 210. 
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Referring again to Fig. 3A, the next step is to select a search area 215 
of reference frame 210 (Step 304). • Search area 215 of Fig. 2 is defined by a 
search range, ±S, and is based on the location of current block 240 in current 
frame 220. Search range S may be predetermined or may be determined 
dynamically. Search range S may be the same for all frames or may depend 
on the time difference between the current frame and the reference frame. 
Search range S can be any positive integer value greater than zero, but less 
than the frame size. Although larger search sizes generally produce a more 
accurate motion vector, a search range S less than the full frame size is 
recommended because the operational complexity of the search may grow in 
proportion to the search area. For the purposes of illustration only, the 
following example assumes the search range is defined as 8 pixels, or S= 8. 
Fig. 2 shows search area 215 defined by S = 8; the point marked (x. y) refers 
to the upper left pixel of current block 240. Assuming the search range 
equals +8. the search area for point (x, y) is the area bounded by (x-8, y-8). 
(x+8. y-8), (x+8, y+8), and (x-8. y+8) determined by the area S pixels from ... 
starting point {x,y) in reference frame 210. The search range in the x and y ■ 
directions can be different, however, for ease of explanation the example 
herein assumes the search range is the same in both x and y directions. 

In a conventional one-pixel-step full search, current block 240 would be 
compared with each possible 8x8 candidate block in the search area. For 
example, current block 240 would first be compared to the 8x8 block in the 
upper leftmost corner of search area 215 and the difference would be 
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calculated. After the first comparison was computed, current block 240 would 
be compared with the 8x8 block beginning one pixel to the right, and the 
difference would be calculated. Following the comparison of current block 
240 with the candidate block beginning at (x+8, y-8), the process would 
continue on the next "row", comparing current block 240 with the candidate 
block with upper left corner (x-8, y-7). Using sum of absolute distortions 
(SAD) as the matching criteria, for example, the full search process may be 
defined as follows: 

7 7 

jSAD{i.J) = J2E\ r(x.y)-s(x^i,y*J) \ 

x=Q y-0 

for -8 <= i <= 8 and -8 <= j <= +8, where r is the current block and s is the 
candidate block. The (ij) pair that produces a minimum value for SAD (i,j) 
defines the motion vector of cun-ent block 240. A motion vector of (1 ,1 ), for 
example, means that a block in the reference frame one pixel horizontally to 
the right and one pixel vertically below the corresponding location of cun-ent ^ 
block 240 in the reference frame closely resembles current block 240. The 

array of values, r(x . y) - s(x + 1, y +j) for x=0, 1 7 and y=0. 1 . :.. 7 are 

recorded as the residue for that current block and candidate block. 

The present invention reduces the number of block comparisons by 
selectively choosing to calculate only a subset of the comparisons calculated 
by the full search. Thus, the process performs a search in one of the 
horizontal or vertical directions. Then, the process performs a limited search 
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in the other direction based on the results of the search in the first direction. 
The order of directions in which the first two searches are performed is 
immaterial, that is, the search may be performed first in the vertical direction 
then in the horizontal direction or vice versa. For the purpose of illustration 
only, the following example assumes that the process begins with a horizontal 
search in the x direction. 

Consistent with the present invention, search area 215 is divided into a 
plurality of subsets of rows (Step 306 of Fig. 3A) as shown in Fig. 4. Each 
subset, or "strip", is Z rows of pixels wide. Z can be any integer greater than 
one and less than or equal to S, the search range. Preferably, however, Z is 
an integer factor of S. In the embodiment depicted in Fig. 4, Z = 8. Search 
area 215 is therefore partitioned into (S/Z) x 2 + 1 = (8/8) x 2 + 1 = 3 strips. 

The process performs the x-directional search by comparing current 
block 240 to a subset of candidate blocks in each strip. In each comparison, 
the y-coordinate of the candidate block will be an integer multiple of Z, or in 
this case, 8. The process calculates the difference from current block 240 : 
and each candidate block using any suitable error calculating Or correlation 
method (Step 308). 

If the difference for a candidate block is less than the difference 
calculated so far for other candidate blocks in this strip or if this is the first 
candidate in the strip, then this candidate is the "best" candidate block (Step 
312). The coordinates and difference of the best candidate block are saved 
for future comparisons (Step 310). If the difference for this candidate block is 
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greaterthan previous difference calculations, the candidate block coordinates 
are discarded. 

If the candidate block is not the last block in the strip (Step 316), the 
process shifts D prxels in the x direction (step 314). Current block 240 is 
compared to the candidate block in the new position (Step 308) and the 
process continues from step 308. D can be any integer value ranging from 1 
to S. If D=1 , for example, current block 240 is compared to every candidate 
block in the strip. 

At the end of a strip (Step 316). the process saves the candidate block 
for that strip by, for example, recording the x and y coordinates that 
correspond to the upper left corner of the candidate block that resulted in the 
least difference from current block 240 (Step 31 8). The process is repeated 
for each strip in search area 215. At the conclusion of the x-directional 
search step (320), the process has recorded one "best" candidate block for 
each strip. The term "best" refers to the candidate block that results in the 
least difference from current block 240. 

Next, the process searches in the other orthogonal direction which, in 
this case, is the y-direction. For the best candidate block of each strip, the 
process defines limited search areas (Step 322, Fig. 3B). Fig. 5 shows 
search area 215 and best candidate blocks. 561 , 562, and 563. Each limited 
search area comprises the partial column formed by the candidate block and 
including some number of additional rows in each direction. In other words, 
each limited search area in the second step is the area formed by sliding the 
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best candidate blocks obtained in step one some number of rows in either 
direction. The number of rows in each direction may be different, however, 
for the sake of explanation, each limited search area includes Y rows in each 
direction. Y may be any integer value greater than one, however, if Y^Z, the 
process will unnecessarily duplicate some difference calculations performed 
in the x-direction. For the purposes of example only, assume Y=Z/2, which in 
this example means Y=4. 

Referring again to Fig. 3B, current block 240 is compared to a first 
candidate block (Step 330). To make the comparison, the left edge of current 
block 240 is aligned' with the left edge of the candidate column. The process 
calculates the difference from current block 240 and the candidate block 
using any suitable error calculating or correlation method. 

If this is the first difference calculated or the difference for this 
candidate block is less than the difference calculated so far for other 
candidate blocks, the coordinates and the difference of the candidate block 
are saved for future comparisons (Step 334). If the difference for this : :: ;};. 
candidate block is greater than previous difference calculations, the candidate 
block coordinates are discarded and the process continues with step 336. 

If the candidate block is not the last block in the limited search area 
(step 336). the process shifts D2 pixels in the y direction (step 335). Current 
block 240 is compared to the candidate block in the new position (Step 330). 
The process continues from step 330. D2 can be any integer value ranging 
from 1 to Y. 
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Referring again to Fig. 3B, if the process reaches the end of a limited 
search area (Step 338). the process goes to the next limited search area and 
continues from step 330. At the end of the limited search areas, the process 
determines the motion vector using the coordinates of the "best" candidate 
block over all the limited search areas (Step 350). The array of pixel 
differences is recorded as the residue for the current block and the "best" 
candidate block. 

B. Sparse Search using integral projection data 

Conventional motion estimation techniques such as block-matching 
algorithms use spatial domain pixel information for estimating the motion 
vectors. When comparing an NxN current block to an NxN candidate block 
using pixel information and a difference calculation such as SAD described 
above, NxN calculations are required. It is possible, however, to compare; • 
the two blocks using other information, such as integral projections. An V 
integral projection of pixel information is a sum of some numberof image S 
pixel values along a certain horizontal or vertical direction." Fig^SA shows a^^^^^^^ 
vertical integral projection; Fig. 6B shows a horizontal integral projection. 

An alternate embodiment of the present invention uses integral 
projection information in conjunction with the sparse search described above 
to further reduce the number of computations performed in determining 
motion vectors. Integral projection information can be obtained by. calculating 
one-dimensional integral projection arrays from pixel information or by 
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calculating integral projection information from discrete cosine transform 
(DCT) coefficient data from, for example, a digital video decoder. In 
particular, the present invention may be used to reduce the computations 
necessary, and therefore decrease the computation time needed, to convert 
from one digital video format to another such as, for example, from digital 
video (DV) format to MPEG format. 

The DCT may be described mathematically as follows: 



Qih.v) = C(h) C(v) ^ P(x.y) COS ( IlIlS^x^] ^os f nXl^yiH] 
yo r.o \ ^6 ) ^ 16 j 



for h=0,1,...7 and v=0.1,...7 and where 



C(/i) = 



2v/2 



/> = 0 



-. />=1 to 7 
2 



and C(v) = 



V = 0 



-, v=1 to 7 
2 



P (x,y) is an 8x8 block of pixel data and Q(h,v) is an 8x8 block of DCT 
coefficients. Mathematically, if an 8 x 8 block of DCT coefficients is used as 
input to a one-dimensional horizontal inverse discrete cosine transform 
(IDCT), the result will be an 8x8 array, the first row of which contains the 
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sums of each column of pixel information multiplied by a constant. The first 
row therefore equals a one-dimensional vertical projection of that 8x8 block of 
pixel information. Similarly, if the same 8x8 block of DCT coefficients is input 
to a one-dimensional vertical I DCT. the first column of the result will be 
equivalent to a one-dimensional horizontal projection of the underlying 8x8 
block of pixel information. The vertical and horizontal integral projections may 
be used to estimate the motion vectors. 

A method consistent with the present invention using integral 
projections is similar to the embodiment described above. Referring again to 
Fig. 3A, the method 'first selects cun-ent block 240 of current frame 220 as 
shown in Fig. 2 (Step 302). Current frame 220 is divided into a plurality of 
identical size search blocks, each comprising NxM pixels. . N and M mayy ^^ 
represent any integer values greater than one, however, when using integral 
projection data obtained by performing a one-dimensional IDCT on the DCTiv 
coefficient array, it is recommended to choose N and M equal to integer^^^^^^ 
multiples of the size of the blocks used during the DCT phase of encoding>Si^ 
For example, if the video encoder performed DCT compression on 8x8 blocks 
of pixel information, choosing N=M=8 is recommended. - 

Referring again to Fig. 3A, the next step is to select a search area 215 
of reference frame 210 (Step 304). The search range. S, may be any integer 
value, however, when integral projections are obtained from 8-point DCT 
coefficients from the video decoding process, it is preferable to select a 
search range that is an integer multiple of 8. The following example assumes 
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the search range is defined as 8 pixels, or S= 8. Assuming the search range. 
S, equals ±8, the search area |s the-area bounded by (x-8, y-8), (x+8. y-8), 
(x+8. y+8). and (x-8, y+8). Search area 215 is then divided into strips as 
shown in Fig. 4, each strip comprising some number, Z, of rows. In this 
embodiment, each "strip" is 8 rows wide so as to be able to take advantage of 
the integral projection data obtained from DCT coefficients. Search area 215 
is therefore partitioned into (S/Z) x 2 + 1 = (8/8) x 2 + 1 = 3 strips. 

The process perfomns the x-directional search by comparing current 
block 240 to a subset of candidate blocks in each strip, as shown in Fig. 4. 
In each comparison; the y-coordinate of the candidate block will be an integer 
multiple of Z. or in this case, 8. The process calculates the difference from 
cun-ent block 240 and each candidate block using any suitable error 
calculating or correlation method (Step 308). 

If the difference for a candidate block is less than the difference 7 ' 
calculated so far for other candidate blocks In this strip or this is the first MP^- 
candidate in the strip, then this candidate is the "best^candidate block (Step::; 
312). The coordinates and difference of the best candidate block are saved 
for future comparisons (Step 310). The process calculates the difference 
from current block 240 and the candidate block using any suitable error 
calculating or correlation method (Step 308). In the first embodiment, the 
calculations were performed on pixel data. In an embodiment of the present 
invention using integral projection data, the one-dimensional integral 
projection array is used instead. If the x-directional search is performed using 
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vertical projections, as the candidate block moves from left to right within a 
strip, most vertical projections for each candidate block can be obtained from 
the block to the left without recalculation. 

If integral projection information is being calculated from pixel data, the 
5 process calculates, N vertical projections for the NxN candidate block and N 

vertical projections for the NxN current block 240. The resulting sums may 
then be used in the difference calculation. For example, using sum of 
absolute distortions (SAD) as the matching criteria, the vertical projection 
: n difference is then calculated as follows: 

M S/\D^(/.y)=X;| r(x)-s(x+;.y) I 

\Q " '0 . . . , . 

J ■'0 Where -S </<+S,-S<y<+S. and; Is an integer multiple of Z. In addition, 

Q I'M 's the vertical projection for the xth column of the current block and s(x+i,j) 

's the vertical projection or sum of the fx+z^th column of the candidate block 
starting at row/ - • ■ 

If the vertical projections for the candidate block are obte'lned from 
"•^ DCT coefficients, vertical projection s(x+i.j) may be computed from two 

blocks. DCT coefficient arrays are calculated for 8x8 blocks such as Block #1 
and Block #2 in Fig. 7. If the coordinates of a candidate block coincide 
exactly with an 8x8 DCT coefficient array, the vertical projection arrays can be 
easily obtained by performing a one-dimensional IDCT on row 0 of the 8x8 
20 DCT coefficient array. If a candidate block straddles two 8x8 DCT arrays 
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horizontally as shown in Fig. 7, vertical projection information is obtained by 
calculating a one-dimensional IDCT for row 0 of the two neighboring 8x8 
arrays for which there are DCT coefficients. 

If a candidate block straddles two 8x8 DCT arrays vertically as shown 
in Fig. 8, the total difference in vertical projection is obtained by calculating 
the weighted sum of the difference values for the two neighboring 8x8 arrays 
for which there are DCT coefficients. For example, the difference value of 
candidate block 820 with upper leftmost point equal to (0,5) as shown in Fig. 
8 is equal to a weighted sum of the difference between block #1 and current 
block 240 and the difference between block #2 and current block 240. 
Candidate block 820 overlaps 3/8ths of block 1 and 5/8ths of block 2, 
therefore, the SAD for candidate block 820 is as follows: 

SAD^iO.5) =3/8 I r(x) -s(x.O) | + 5/8 J) j r(x) -s(x;8) | ■ . > 




where r^x|ls the vei^ical projection for the xth column of the current block, ^^^^^^ 
s(x,0) is the vertical projection or sum of the xth column of block #1 and s(x,8) 
is vertical projection or sum of the xth column of block #2. ■ 

Referring again to Fig. 3A, if this is the first difference calculated or the 
difference for this candidate block is less than the difference calculated so far 
for other candidate blocks in this strip (step 312), the coordinates and 
difference of the candidate block are saved for future comparisons (Step 
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310). If the difference for this candidate block is greater than previous 
difference calculations, the candidate block is discarded. 

If the candidate block is not the last block in the strip (Step 316), the 
process shifts D pixels in the x direction (step 314). Current block 240 is then 
compared with the candidate block in the new position (step 308) and the 
process continues from step 308. D can be any integer value ranging from 1 
to S. If D=1 . for example, current block 240 is compared to every candidate 
block in the strip. 

At the end of a strip, the process records the x and y coordinates that 
correspond to the upper left comer of the candidate block that resulted in the 
least difference from current block 240 (Step 318). The process is repeated 
for each strip in search area 215. . At the conclusion of the x-directional:^ ; 
search step, the process has recorded one "best" candidate block for each 
strip (Step 320). The temi "best" refers to the candidate block that results i 
the least difference from cun-ent block 240. ; /5 ' ' - f ^^'/-^tr^y^^tip'^-^ 

Next, the process searches In the y-direction;4Fig. 5 shows searchl^ ^^ 
area 215 in reference frame 210 and best candidate blocks 561 ,^ 562. an3 : 
563 from the x-directional search. For the best candidate block of each strip, 
the secondary limited search area comprises the column formed by the best 
candidate block and includes Y rows in each direction. Y may be any integer 
value greater than one, however, if Y^S, the process will unnecessarily 
duplicate some difference calculations performed in the x-direction. For the 
purposes of example only, assume Y=S/2 or in this case, Y=4, 
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The process can shift any number of pixels (1 to Y) in the y direction, 
however, for ease of explanation, we assume the process shifts one pixel in 
the y direction each time. The compare process may begin, for example, by 
comparing current block 240 and the candidate block with upper leftmost 
point equal to (m. n-Y) if the upper leftmost point of best candidate block 561 
is (m. n) (Step 330). 

The process calculates the difference from current block 240 and the 
candidate block using any suitable error calculating or correlation method 
(Step 330). In the first embodiment, the calculations were performed on pixel 
data. In an embodiment of the present invention using integral projection 
data, the one-dimensional integral projection array is used instead. 
Performing the second search process, the y-directional search in this . . ■ :/ 
example, using both horizontal and vertical projections may produce more 
accurate motion estimation. Using horizontal integral projections only, 
^ however, reduces the number of operations performed and may be preferable 
in some applications, -v:- •v'^ ; • > - ; 'f :^rv^'?^^':^'\f'^M^.\ 

integral projection information is being calculated from pixel data, the 
process calculates N horizontal projections for the NxN candidate block and N 
horizontal projections for the NxN current block 240. The resulting sums may 
then be used in the difference calculation. For example, using sum of 
absolute distortion (SAD) as the matching criteria, the horizontal projection 
difference is then calculated as follows: 




where n-Y < / < n+Y, where t(y) is the horizontal projection for the yth row of 
the current block and u(m. j+y) is the horizontal projection or sum of the yth 
row of the candidate block. 

The vertical projection difference can be calculated using: 

7 

jSAD^im. j)=J2\ r(x) -s(x + m, y) | 

Jf=0 

Where -S <j < + S. where r(x) is the vertical projection for the xth column of 
the current block and s(x + m,j) is the vertical projection of the xth column of 
the candidate block. The total difference is the sum of the horizontar' . . 4 , V 
projection difference and the vertical projection difference. W 

: ; ; I projection data for the candidate blocks are obtained from DCT ' 
coefficients, the total difference calculation Is more complicated;^!ln th^x-l^^^^^ 
directional search, if the upper edge of current block 240 is aligned with the 
upper edge of the candidate row, the y-coordinate will be an integer multiple 
of Z, or in this case, 8, making easy the calculation of integral projections from 
DCT coefficients. In the second step of the search, whether vertical or 
horizontal, the candidate blocks will not be aligned with the 8x8 arrays of DCT 
coefficients. For example, the total difference of candidate block 910 with 
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upper leftmost point equal to (5,2) as shown in Fig. 9 is equal to a weighted 
sum of the differences between blocks #1. #2, #3, #4. and current block 240. 
For example, the weighted sum of the differences of the four neighboring 
blocks can be calculated as ((8-2)/8)*(the difference from search using 
vertical projections at (5,0)) + (2/8)*(the difference from search using vertical 
projections at (5,8)) +((8-5)/8)*(the difference from search using horizontal 
projections at (0,2)) + (5/8)*(the difference from search using horizontal 
projections at (8,2)). 

Mathematically, the SAD for candidate block 910 is as follows: 

J 

S/\D(5.2)=i^J^| r(x)-s(x+5.0) j + - ^ j r(x) -s(x+5, 8) | 

^^~l:\ny)-u{0.y*2) \^^j2\ t(y)-u(B,y^2) \ 
o y=o 8 y=o 

where r(x) is the vertical projection for the xth column of the current block, t(y) 
is the horizontal projection for the yth row of the current block, s(x+5,0) and 
s(x+5,8) are the vertical projections or sums of the (x+5)ih column of block #1 
and block #3. respectively, for 0 s x :s 2; s(x+5,0) and s(x+5, 8) are the 
vertical projections or sums of the (x+5-8)ih column of block #2 and block #4, 
respectively, for 3 s x s 7; u(0,y +2) and u(8, y+2) are the horizontal 
projections or sums of the (y+2)ih row of block #1 and block #2, respectively. 
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for 0 5 y s 5; u(0, y+2) and u(8. are the horizontal projections or sums of 
the fy+2-8;th row of block #3 and block,#4, respectively, for 6 s y ^ 7. 

The process increases the vertical index j with each comparison until it 
finds the best j that produces the minimum difference between the candidate 
block and current block 240. This procedure is repeated for each of the mini 
search areas defined by best candidate blocks, 561. 562. and 563 and a best 
candidate block is recorded. With each comparison, if the candidate results 
in the minimum total difference from current block 240. the coordinates and 
the total difference of the candidate are recorded. When all the mini searches 
are complete, the recorded best candidate block is the reference block that 
best resembles current block 240. Finally, the process calculates the motion 
vector and residue from the best candidate block (m.l) and current block 240. 

Calculating integral projections using 1-D IDCTs of frequency domain 
information requires approximately 1 1 multiplications and 29 additions for 
each set of 8 projections. If the projections are calculated purely in the 
spatial domain, it takes about 56 operations for each set of 8 projections. The 
performance gain of the frequency domain method will be higher if only the 4- 
point IDCT of the four lower frequency coefficients is used. This process is 
equivalent to performing further filtering and decimation on the 8 projections. 
In this case, only four integral projections in every 8x8 block are used. The 
four projection data obtained from the 4-point IDCT are equivalent to the 
sums of every two adjacent integral projections as illustrated in Fig. 10A. It is 
also possible to use 2 projections in every 8x8 block by using 2-point IDCT on 
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two lower frequency coefficients or performing another step of summation 

similar to the procedure in Fig.. 1 0B.- Using fewer projections can further 

reduce the number of operations in motion search. 

C. System 

Fig. 11 illustrates a system 1105 consistent with the present invention. 
As shown in Fig. 1 1 , a processor 1 1 10 is connected to at least one 
input/output (I/O) device 1 120 via any suitable data connection. I/O device 

1 120 can be any device capable of passing information to or receiving data 
from processor 1110. By way of example only, I/O device 1 120 may be a 
digital camcoder connected through an IEEE 1394 interface. Processor 1110 
may be any commonly available digital processor such as, for example, a 
Pentium processor. Processor 1110 may be a single processor or multiple 
processors. Faster processors, however, will decrease execution time of the 
invention. 

The system of the present Invention also includes memory 1130 
capable of storing data processed by processor 1110 and data sent to or ; : : 
received from I/O device 1 120. System 1105 may be connected to a display 
1 140, such as a cathode ray tube (CRT), for displaying information. 
Processor 1110, I/O device 1120, memory 1130, and display 1140 are 
connected via a standard system bus 1160. Fig. 11 shows an exemplary 
network where each hardware component may be implemented by 
conventional, commercially available computer systems components. 
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Fig. 12 illustrates a processor consistent with the present invention. 
Processor 1110 may comprise.one or nnore memory management units 
(MMU) 1210, one or more processor element arrays 1220, and one or more 
accumulator units 1230. Processor element array 1220 may comprise an 
array of processor elements 1225. Processor elements 1225 may comprise, 
for example, a subtraction and adder units for calculating the SAD between 
the blocks. MMU 1210 may be used to buffer the data for processor element 
array 1220. Accumulator unit 1230 may be, for example, an adder unit that 
adds the outputs from processor elements 1225. 

Consistent with one implementation of the invention, processor 1110 
executes one or more sequences of one or more instructions contained in 
memory 1 130. Such instructions may be read into memory 1 130 from a 
computer-readable medium via Input/output device 1120. Execution of the 
sequences of instructions contained in memory 1130 causes processor 1110 
to perfonn the process steps described herein. In an alternative 
implementation, hard-wired circuitry may be used in place of or in ' 
combination with software instructions to implement the invention'. Thus 
implementations of the invention are not limited to any specific combination of 
hardware circuitry and software. 

The term "computer-readable medium" as used herein refers to any 
media that participates in providing instmctions to processor 1 1 10 for 
execution. Such a medium may take many forms, including but not limited to, 
non-volatile media, volatile media, and transmission media. Non-volatile 
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media includes, for example, optical or magnetic disks. Volatile media 
includes dynamic memory, such as memory 1 130. Transmission media 
includes coaxial cables, copper wire, and fiber optics, including the wires that 
comprise system bus 1 160. Transmission media can also take the form of 
acoustic or light waves, such as those generated during radio-wave and 
infra-red data communications. 

Common forms of computer-readable media include, for example, a 
floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic 
medium, a CD-ROM, any other optical medium, punch cards, papertape, any 
other physical medium with patterns of holes, a RAM, PROM, and EPROM, a 
FLASH-EPROM, any other memory chip or cartridge, a carrier wave as . 
described hereinafter, or any other medium from which a computer.can read. • 
Network signals carrying digital data, and possibly program code, to and from 
system 1105 through system bus 1160 are exemplary forms of carrier waves .'- 
transporting the information. In accordance with the present invention, . # 
program code received by system 1 105 may be executed by processor 1 1 10 ' 
as it is received, and/or stored in memory 1 1 30, or other non-voiatile storage 
for later execution. . 

It will be apparent to those skilled in the art that various modifications 
and variations can be made in the methods and systems consistent with the 
present invention without departing from the spirit or scope of the invention. 
The true scope of the invention is defined by the following claims. 



